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We report molecular dynamics simulations of a hydrophobic polymer-chain in aqueous solution 
between 260 K and 420 K at pressures of 1 bar, 3000 bar, and 4500 bar. The simulations reveal a 
hydrophobically collapsed state at low pressures and high temperatures. At 3000 bar and about 
260 K and at 4500 bar and about 260 K, however, a transition to a swelled state is observed. The 
transition is driven by a smaller volume and a remarkably strong lower enthalpy of the swelled state, 
indicating a steep positive slope of the corresponding transition line. The swelling is stabilized almost 
completely by the energetically favorable state of water in the polymers hydrophobic first hydration 
shell at low temperatures. Although surprising, this finding is consistent with the observation of a 
positive heat capacity of hydrophobic solvation. Moreover, the slope and location of the observed 
swelling transition for the collapsed hydrophobic chain coincides remarkably well with the cold 
denaturation transition of proteins. 
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I. INTRODUCTION 

Hydrophobic effects have been shown to be of rele- 
vance for a wide range of physicochemical and biophys- 
ical phenomena 2, 3]. In particular, they are seen as 
an important driving force with regard to the folding of 
proteins Consequently, a wealth of studies of hy- 

drophobic interactions using molecular simulation tech- 
niques have been undertaken over the past three decades 
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I22L l23l 124 l25l l26l l27| . Simulation studies have revealed 
that the contact state of a pair of hydrophobic parti- 
cles in aqueous solution is entropically stabilized at am- 
bient conditions [Tol fllj . In addition, the contribution of 
the solvation heat capacity has been recognized recently 
|2Q, 21, 22, 23,, 24, .25, .26, 27}. In general, the disso- 
lution of hydrophobic particles is accom pan ied with an 
increase of the associated heat capacity |2^ |2^ . As 
a consequence, the dissolution of a hydrophobic particle 
is found to be increasin gly enthalpically stabilized with 
decreasing temperature |23 . 

Biopolymers such as proteins remain stable and func- 
tional only in a limited pressure and temperature range 
[311 133 . Issf . Increasing temperature lead to structural 
changes differing from the native folded state. This is 
often accompanied with large fluctuations and aggrega- 
tion phenomena. Hence, pressure effects on proteins are 
of interest in biotechnology and biology [3j, as pres- 
sure is shifting the equilibrium of protein configurations 
without increasing thermal fluctuations 35, 36, 37, 38]. 
Proteins undergo unfolding upon addition of pressures 
above 2 kbars. High pressures are also able to dissociate 
protein complexes The solvent water plays a 

crucial role in the effect of pressure in protein unfolding 
[33L I36L im li^ and the addition of co-solvents is found 
to have a significant influence on the size, location and 



'Electronic address: "diet max.paschekQudo.edul 



shape of the stability region of proteins |43, |4j| . 

At high pressures (> 2 kbars), the volume of proteins 
upon unfolding decreases. This seemed to be inconsistent 
with the assumption that protein unfolding is equivalent 
to the transfer of hydrophobic groups from the protein 
interior to the aqueous solvent, since the volume change 
upon transfer of hydrophobic groups to water are posi- 
tive. Hummer et al. [45j suggested a scenario in which 
pressure unfolding of proteins is modeled as the trans- 
fer of water into the protein hydrophobic core with in- 
creasing pressure. The transfer of water molecules into 
the protein interior is essential for the pressure unfolding 
process, leading to the dissociation of close hydropho- 
bic contacts and subsequent swelling of the hydropho- 
bic protein interior through insertions of water molecules 
|45| . The characteristic features of water-mediated inter- 
actions between hydrophobic solutes in water are found 
to be pressure-dependent. In particular, with increas- 
ing pressure the solvent-separated configurations in the 
solute-solute potential of mean force is stabilized with re- 
spect to the contact configurations. In addition, the de- 
solvation barrier increases monotonically with respect to 
both contact and solvent-separated configurations. The 
locations of the minima and the barrier move toward 
shorter separations, and pressure effects are considerably 
amplified for larger hydrophobic solutes 0, 0, 0| . 

Pressure also changes the entropy/enthalpy balance of 
the hydrophobic interactions. Ghosh et al. found that 
the contact minimum is dominated by entropy, whereas 
the solvent-separated minimum is stabilized by favorable 
enthalpy of association 0, Ts^ . Both the entropy and 
enthalpy at the contact minimum seem to change lit- 
tle with increasing pressure leading to the relative pres- 
sure insensitivity of the contact minimum configurations. 
In contrast, the solvent-separated configurations are in- 
creasingly stabilized at higher pressures by enthalpic con- 
tributions that prevail over the slightly unfavorable en- 
tropic contributions to the free energy [l^ . 

In this contribution, we focus particularly on the sce- 
nario proposed by Hummer et al. lisl of water penetrat- 



2 



ing into the protein interior at elevated pressures. We 
study, however, a very much simpHfied model system of 
protein: A fully hydrated polymer-chain, consisting of 
20 interconnected hydrophobic particles. The polymer- 
chain approach has been recently advocated by Chandler 
and co-workers [ifilll^li^ . suggesting that the collapse of 
a hydrophobic polymer chain is driven by a drying tran- 
sition. In Ref. 47] the hydrophobic chain is modeled by 
repulsive polymer/water interactions only, whereas the 
solvent is represented by a coarse grained model. In the 
recent work of Ghosh et al. 49j , molecular dynamics sim- 
ulations of a polymer-chain in an explicit solvent reveal 
the effect of salt concentrations on the polymer configura- 
tion, increasingly favoring compact folded configurations 
of the polymer with increasing salt concentration. In 
the above mentioned studies a "stiff" polymer chain was 
employed, which exhibits a stretched equilibrium config- 
uration in absence of a solvent. Thus the energy needed 
to bend the polymer-chain is used to counterbalance the 
tendency to minimize the solvent accessible surface. In 
the present study no such terms are employed, hence a 
chain of linked hydrophobic particles, free of bond-angles 
and dihedral potential barriers is considered. Thus the 
reported configurational changes of the polymer have to 
be completely attributed to the influence of the solvent. 
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II. COMPUTATIONAL METHODS 



MD Simulation details 



TABLE I: Simulation protocol for the performed MD-simulations. 
t: Simulation time; (V): Average box- volume; {E): Average total 
energy. 



We report molecular dynamics (MD) simulations of a 
purely hydrophobic polymer-chain dissolved in an aque- 
ous solution. The polymer-chain consists of 20 hydropho- 
bic polymer beads, represented by Lennard-Jones inter- 
action sites with axx=3.975A, exx fc^^=214.7K [26l| . 
The water-polymer cross parameters were obtained us- 
ing the conventional Lorentz-Berthelot mixing rules with 
) /2 and e^j = y/euejj. The polymer sites 
are linked by rigid bonds of 4.2 A length. This bond- 
length was determined to represent the bead-bead con- 
tact distance for non-linked particles in aqueous solution 
pyf. In addition, this value corresponds roughly to the 
distance between adjacent hydrophobic Valine sidechains 
in a polypeptide. All intramolecular non-bonded in- 
teractions were taken into account, except interactions 
between adjacent bonded sites. No additional bond- 
bending or torsional potentials were used. The water 
phase is represented by 1000 TIP5P water molecules 
|50|. The TIP5P model was chosen, since it represents 
the hydrophobic solvation behavior of water on the low- 
temperature side probably most realistically among the 
simple point charge models ^26J. Moreover, the temper- 
ature dependent strength of the hydrophobic interaction 
was found to be quite critically linked to the tempera- 
ture dependence of waters expansivity p^ . The simula- 
tions discussed here were carried over a broad tempera- 
ture range at pressures of 1 bar, 3000 bar, and 4500 bar. 



Individual MD-simulations extend up to 100 ns, while the 
total simulation time adds up to about 0.76 ^s. For com- 
pleteness, a detailed simulation protocol is given in Table 
III 

The MD-simulations are carried out in the NPT en- 
semble using the Nose-Hoover thermostat 51, 52] and 
the Rahman-Parrinello barostat 53, 54] with coupling 
times tt ~ 1.5 ps and — 2.5 ps (assuming an isother- 
mal compressibility of xt = 4.5 10~^bar~^), respec- 
tively. The electrostatic interactions are treated in the 
"full potential" approach by the smooth particle mesh 
Ewald summation [55j with a real space cutoff of 0.9 nm 
and a mesh spacing of approximately 0.12 nm and 4th 
order interpolation. The Ewald convergence factor a 
was set to 3.38 nm~^ (corresponding to a relative accu- 
racy of the Ewald sum of lO"'"^). A 2.0 fs timestep was 
used for all simulations. Solvent constraints were solved 
using the SETTLE procedure [11], while the SHAKE- 
algorithm was used for the polymer constraints |57j ] . For 
all simulations reported here the GROMACS 3.2 pro- 
gram 58, 59] was used. Statistical errors in the analysis 
were c omp uted using the method of Flyvbjerg and Pe- 
tersen |60j. For each system an initial equilibration run 
of about 1 ns length was performed using the Berend- 
sen weak coupling scheme for pressure and temperature 
control (tt = rp = 0.5 ps) |61|. 



3 



It should be mentioned that the completely stretched 
polymer extends to about 8 nm, which exceeds the used 
box length about 2.5 times. However, in practice, con- 
tacts between the polymer and its virtual image were not 
observed during the simulation runs discussed here. 



B. Energy Analysis 

In order to assign potential energies to individual 
molecules and thus to be able to distinguish between 
contributions from the "hydration shell" and from the 
"bulk" , we determine energies by a reaction field method 
based on the minimum image cube. This "cubic" cut- 
off procedure has been originall y p roposed by Neumann 
. Roberts and Schnitker [63Ll64| have shown that the 
obtained energy estimates are comparing very well with 
the Ewald-summation including tin-foil boundary condi- 
tions. 

For convenience, we divide the potential energies in 
contributions assigned to the individual molecules with 

M 

E = Y.E, 

(l ^ \ 
~ I o5Z^*^ + Ei^corr. , (1) 

Vs=i ; 

where Ei is the potential energy assigned to molecule 
i, M is the total number of molecules. The molecule- 
molecule pair energy 




is then obtained as sum over discrete interaction sites 
a and /3, with riajis = \fji3 — fia\ based on the 
molecule/molecule center of mass minimum image sep- 
aration. We employ long range corrections Ei^corr. = 
^ijcorr. + ^ijorr. accountiug for electrostatic, as well as 
Lennard Jones interactions. The electrostatic correction 

is a reaction field term, corresponding to the cubic cut- 
off, assuming an infinitely large dielectric constant. Here 
'^i = Qiafia IS the dipole moment of molecule i, 
D = di is the total dipole moment of all molecules in 
the simulation cell and V is the instantaneous volume of 
the simulation box. Finally, also the long range Lcnnard- 
Jones corrections for to the minimum image were taken 
into account, as outlined in Ref. p^ . 



III. RESULTS 

The structure of the dissolved hydrophobic polymer- 
chain is characterized by its radius of gyration Rq =^ 
1/20 in - c)^ with c = l/20X;^°if,. The typical 
conformation of this polymer at high temperatures and 
low pressures is a compact, collapsed state with a Rg be- 
tween 0.5 nm and 0.6 nm, as shown in Figure|3| A snap- 
shot of a representative coUapsed-chain configuration is 
shown in Figure^ The conformational distribution with 
respect to Rq is found to be narrow, with a half width at 
half maximum of Rg of about 0.1 nm. Test-simulations 
at BOOK and ambient pressure conditions, starting with 
a swelled configuration of Rg ~ 1 nm, show a collapse on 
a timescale < Ins. The polymer/water center of mass 
pair correlation functions (given in Figure 0J) reveal that 
in the collapsed state water is completely excluded from 
the polymer interior. 

The temperature dependence of the average radius of 
gyration at a pressure of 1 bar is characterized by a shal- 
low minimum at about 300 K. With increasing temper- 
ature the distribution of Rg more or less maintains its 
shape, but is becoming broader with its maximum shift- 
ing to larger values as shown in Figure [3 The high- 




FIG. 1: Representative configurations of the hydrophobic 
polymer in aqueous solution as obtained from simulations 
3000 bar. Top: Collapsed configuration observed at 300 K. 
Bottom: Swelled configuration observed at 260 K. 
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Bulk: 260 K; 3000 bar 280 K; 4500 bar 

£(H20)/kJmor^ -43.657 ± 0.004 -42.081 ± 0.003 

V'm(H2 0)/cm^mori 16.02 ± 0.05 15.50 ± 0.02 

Hydration Shell: 260 K; 3000 bar 280 K; 4500 bar 

"collapsed" "swelled" "collapsed" "swelled" 

iV(H20) 126.5 ±2.1 184.2 ±2.8 129.4 ±1.3 173.0 ± 2.0 

V(Shell)/nm^ 4.51 ± 0.09 5.96 ± 0.08 4.49 ± 0.03 5.58 ± 0.07 

V"™(H2 0)/cm^mor^ 21.47 ± 0.05 19.49 ± 0.05 20.97 ± 0.07 19.46 ± 0.07 

£;(H20)/kJmor^ -43.63 ±0.05 -44.43 ± 0.03 -41.85 ±0.02 -42.34 ±0.03 

S(Polymer)/kJmol-i -149.7 ±0.5 -141.3 ±0.6 -141.2 ±0.2 -126.7 ±0.5 



TABLE II: Energies and volumes obtained for the bulk and hydration shell for the given statepoints. The "hydration shell" is defined as 
the volume with a distance i?<0.52nm to any polymer site, whereas the "bulk" is obtained for distances of i?>1.0nm. y(Shell) denotes 
the corresponding volume occupied by the polymer and its first hydration shell. Af(H20) gives the average number of hydration shell 
water molecules. The hydration shell water molar volumes are according to Kn(H20) =V(Shell)/Af (H2O) * Na, where Na is the Avogadro 
number. The potential energies E of the polymer and water in bulk and shell are given per molecule, "collapsed" and "swelled" states 
are defined by Rq <0.fi5nm and Rq >0.8nm, respectively. 
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FIG. 2: Time-evolution and corresponding probability density distribution of the polymers radius of gyration. Shown are the 
lowest temperatures at 1 bar (bottom), 3000 bar and 4500 bar (top). The dashed line indicates the probability distribution 
obtained for the highest temperatures, respectively. 



temperature behavior is more or less similar for all pres- 
sures discussed here. In addition, a slight penetration of 
water into the polymer interior is observed for 420K (see 
Figure 

At lower temperatures, however, a different behavior 



is starting to emerge. From the time evolution of Rq, 
shown in Figure |21 it is evident that at about 260 K and 
1 bar the polymer increasingly starts to explore extended- 
chain configurations. These extended-chain states are oc- 
curring only rather infrequently and are short-lived with 
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a life-time of about 1 ns. We would like to emphasize 
that the 280K/3000bar and 300K/4500bar trajectories 
show a similar behavior. Again, the extended-chain con- 
figurations are rather short-lived with the chain quickly 
returning to the collapsed state. In a quantitative man- 
ner, an increase of the population of swelled states is 
found, which is due to the apparent higher frequency of 
large amplitude Rq -OxLCtuations. The most striking dif- 
ference, however, is observed for 260 K at 3000 bar and 
280 K at 4500 bar . Here the extended-chain configura- 
tions are dominating, although an equilibrium between 
collapsed and swelled states is still maintained. The ra- 
dius of gyration is showing an apparently bimodal distri- 
bution at 260 K at 3000 bar and a broad distribution at 
280 K at 4500 bar. The representation of both (swelled 
and compact) at each of these states, strongly suggests 
that the corresponding conformational transition tem- 
peratures are located in close proximity to the tempera- 
tures indicated. The temperature dependence of Rq at 
different pressures, as given in Figure El suggests that 
the high-temperature coUapsed-chain state is more com- 
pact at elevated pressures, which would correspond to 
a more compressed coiled state. At 280 K and 3000 bar 
and 260 K and 4500 bar, however, the situation has al- 
ready changed, and the tendency to explore more ex- 
tended configurations leads to an increase in Rq which 
is quickly progressing upon cooling. The low tempera- 
ture destabilization of the collapsed state is in line with 
the decreased stability of hydrophobic contacts observed 
for elevated pressures 18, 19, 45]. Figure 13 implies that 
the transition towards a swelled state at 3000 bar and 
4500 bar occurs in a rather narrow temperature interval, 
suggesting a rather large enthalpy difference between col- 
lapsed and swelled state. In other words: it shows a large 
"cooperativity" . Thermodynamical consistency requires 
that the swelled low temperature state has to be energet- 
ically more stable than the compact state and that it has 
to occupy a smaller total volume. Figure[51shows a super- 
position of the time-evolution of the potential energy and 
box- volume, as well as Rg for 260 K at 3000 bar. In order 



to make the trends more clearly visible, the noise in the 
fluctuations have been reduced by Savitzky Golay filter- 
ing [s^. Both, the potential energy as well as the box- 
volume are clearly anti-correlated with respect to Rg- 
The swelled state of the polymer is enthalpically stabi- 
lized with energy difference of AEu = ii^(swell.) — E{co\\.) 
of -182kJmor^ at 3000 bar and of -82 kJmol"^ at 
4500 bar. The swelled state leads also to smaller box vol- 
umes, with AVu = y(swen.) - V{co\\.) of -20 mlmol"^ 
at 3000 bar and of — lOmlmol"^. We would like to point 
out the observed enthalpy and volume changes are of 
similar magnitude as observed for some proteins [3^ and 
peptides |6a| . 

In the following paragraphs we would like to show 
that the energy and volume change can be attributed 
almost quantitatively to the first hydration shell of the 
polymer and have to be almost exclusively attributed to 
the solvent. Figure El depicts the polymer/water-oxygen 
site-site pair correlation functions for the collapsed and 
swelled states obtained at 260 K and 3000 bar. The first 
minimum, located at a distance of 0.52 nm, indicates the 
spatial dimension the first hydration shell. The aver- 
age number of water neighbors around each polymer site 
changes about 40% from 11.5 to 16 upon swelfing (at 
3000 bar) . At 4500 bar a qualitatively similar behavior is 
observed. In order to quantify the changes in the hydra- 
tion shell, we define the volume occupied by the polymer 
and the first hydration shell as the volume with a dis- 
tance of i? < 0.52 nm with respect to any polymer-site. 
The properties obtained for the hydration shell of the 
collapsed and swelled polymer and the water bulk are 
summarized in Table ITll 

The polymer/water center of mass pair correlation 
function as well as the potential energy of water as a 
function distance to the polymer center for the collapsed 
and swelled states are given in Figure From Figure 
[3it is evident that for the region of the polymer/ water 
interface (r«0.9nm) water is 0.3kJmol~^ more stable 
as in the bulk. However, due to the lack of water neigh- 
bors, a further penetration of individual waters into the 
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FIG. 3: Average radius of gyration Rg of the hydrated poly- 
mer as obtained for all temperatures at 1 bar, 3000 bar and 
4500 bar. The lines are drawn to guide the eye. 
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FIG. 4: Polymer/ water center of mass pair correlation func- 
tions for all simulated state points. 
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FIG. 5: Time evolution of the radius of gyration of the poly- 
mer (black) , potential energy (red, top) and volume (red, bot- 
tom) as obtained from the simulation at 260 K and 3000 bar. 
Potential energies and box-volumes were directly taken from 
GROMACS simulation output. In order to reduce the noise, 
all data were smoothed using with the same Savitzky-Golay 
filter EBi- 
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FIG. 6: Polymer-bead/water-oxygen site-site pair correlation 
function 5X-Ow('') ^^'^ integrated number of nearest neigh- 
bors 7V„ (r) for the collapsed and swelled states at 260 K and 
3000 bar. 



collapsed hydrophobic coil is energetically unfavorable, 
as the steep increase in Figure [7|d indicates. When 
drawing the balance over all water molecules in the 
hydration shell, energy gains and penalties almost cancel 
out completely and the average potential energy of a 
water molecule in the hydration shell of the collapsed 
polymer is nearly identical to the value observed for the 
water bulk at 3000 bar and slightly more unstable at 
4500 bar (see Table Hlj) . For the case of the swelled poly- 
mer, the energy penalty is absent and the water in the 
hydration shell gains — 0.78kJmol~^ potential energy 
per water molecule with respect to the bulk on average 
(— 0.26 kJ mol~^ at 4500 bar). The potential energy 
difference between collapsed and swelled state is about 
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FIG. 7: (a) Polymer/water center of mass pair correlation 
function for the collapsed and swelled states at 260 K and 
3000 bar. (b) Potential energy of the water molecules as a 
function of distance to the polymer center. 
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FIG. 8: Free volume fraction accessible to a small hard 
sphere particle as a function of distance to the polymer cen- 
ter. Given are dependencies for the collapsed and swelled 
states at 260 K and 3000 bar. The diameter of the particle 
with oy-ow ~ 2.5 A and cry_x ~ 2.93 A is scaling as the 
corresponding Lennard- Jones sigmas. 



AK « [i;(shell,swell.) - £;(bulk)] Ar(shell,swell.) - 
[£;(shell,coll.) - £;(bulk)] A^(shell,coU.) + 
i?(polymer, swell.) — i?(polymer, colli.) = — 138kJmoP^ 
(60 kJ mol""'^) which already accounts largely for the 
observed total energy difference of about — 182kJmol~^ 
(-82kJmor^). The potential energy of the polymer 
changes just by 8.4kJmol~^ (14.5kJmol~^), which is 
due to the fact that the loss of intramolecular interac- 
tions when going from the collapsed to the swelled state 
is almost completely compensated by polymer/solvent 
interactions (The values for 4500 bar are given in paren- 
theses). Hence the the extended-chain configurations at 
low temperatures appear to be largely solvent-stabilized. 

The volume change upon swelling is Al^ k, 
[V;„(shen,swell.) - K„(bulk)] 7V(shell,swen.) ' - 
[V;„(shen,con.) - T/(bulk)] 7V(shell,con.) = -57mlmor^ 
at 4500 bar and AKi,„ « -21 mlmol"^ for the 4500 bar- 
isobar, origin of the negative volume change upon 
swelling is indicated in Figure 01 Here the free volume 
fraction accessible to a small sphere is shown as a 
function of distance to the polymer center. For the 
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collapsed state, apparently an excess free volume in the 
polymer interior is available to a hard sphere particle 
which is absent in the swelled state. The hydration 
shell around the hydrophobic chain appears to be more 
tightly packed, so that the average free free volume 
fraction is even lower than for the water bulk. Hence 
the increase in solvent accessible surface (increasing 
number of hydration shell waters) is overcompensated 
by the decrease in molar volume of water in the swelled 
hydration shell state. The hydration shell volume of 
the collapsed state includes, of course, the volume of 
the hydrophobic core, which is made accessible to the 
solvent upon unfolding of the chain. 

Finally, we would like to compare the behavior of the 
hydrophobic polymer to the experimentally obtained sta- 
bility diagrams of proteins. Interpolating the i?(5-data 
from Figure 121 and assuming that the transition appears 
at Rg « 0.7 nm, we obtain ~ 268 K and 282 K for 
3000 bar and 4500 bar, respectively. Using the enthalpies 
AiJ„ « AKtt -I- Pl^Vu and volumes AT4 according the 
data from Table ITU we obtain a slope of the coexistence 
line of adpeq/dT = AHu/T^AVu of about 100 barK"^ 
for both pressures, seemingly consistent with the approx- 
imate transition temperatures. Although having more 
data would be desirable, we might conjecture that the 
slope does not seem to change much, but the enthalpy 
and volume differences tend to decrease with increasing 
pressure. Since the differences might disappear com- 
pletely at higher pressures, the swelling transition of 
the hydrophobic polymer probably ends up in a "critical 
point" (The term "critical point" might not be fully ap- 
propriate here since swelling transition has a monotonous 
character due the finite size of the polymer). In Figure 
l^lthe swelling transition line, as well as the stability di- 
agrams for Staphyloccal Nuclease ^^'^ ubiquitin 
are shown. We would like to point out that the loca- 
tion and slope of the swelling transition shows remark- 
able similarity to the given cold-denaturation lines. The 
hydrophobic polymer seemingly behaves as suggested by 
the water penetration scenario according Hummer et al. 
|45| . In the present case the energy stabilization of the 
swelled configuration is dominated the energy gain of the 
hydration water. In a real polypeptide this is not neces- 
sarily true. Backbone hydration and the equilibrium be- 
tween intra- and intermolecular hydrogen bonds will cer- 
tainly play an important role. How this delicate balance 
might influence the equilibrium between swelled and col- 
lapsed configuration should be further investigated. The 
high temperature side of the protein stability diagram 
might not be accessible by our simple model since it is 
largely related to internal secondary structural transfor- 
mation to a "molten globule" state ||67|. In addition, we 
would like to mention that the observed cold swelling of 
the hydrophobic polymer has some similarity to the sce- 
nario proposed for the swelling of polymeric tropo-elastin 
UTX)n cooling, as advocated by the work of D.W. Urry 

In the present simulation study the energy gain upon 
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FIG. 9: Sketch of the swelling line of the hydrophobic poly- 
mer used in this study compared with experimental stability 
diagrams as obtained for SNase 41] and Ubiquitin |43| . 



unfolding is apparently a consequence of the increasing 
water-water pair interactions in the hydration shell at 
low temperatures |27j . The hydrogen bond network in 
the hydration shell starts disintegrating upon heating 
more strongly as in the bulk since no hydrogen bonds 
can be formed to the hydrophobic Lennard-Jones parti- 
cles [23. The counterbalance between strengthened hy- 
drogen bonds and enhanced disintegration of the hydro- 
gen bond network is widely regarded as the mechanism 
leading to the positive solvation heat capacity associated 
with hydrophobic hydration a consequence, 

the observed transition temperature might be critically 
influenced by the change of the heat capacity of the sol- 
vent in the hydration shell compared to the bulk. This 
might also shed light on the mechanism how co-solvents 
affect the (cold) unfolding transition of proteins. 



IV. CONCLUSION 

Molecular dynamics simulations of a hydrophobic 
polymer-chain in aqueous solution between 260 K and 
420 K at pressures of 1 bar, 3000 bar, and 4500 bar reveal 
a hydrophobically collapsed state at low pressures and 
high temperatures. At about 268 K and 3000 bar and at 
282 K and 4500 bar a transition to a swelled state is ob- 
served. The transition is driven by a smaller volume and 
a remarkably strong lower enthalpy of the swelled state. 
The volume effect is basically due to a smaller net- volume 
of the extended hydrated state compared to the collapsed 
state exhibiting hydrophobic cavity volumes and pene- 
tration of the internal volume by water. Moreover, the 
extended-chain structure is almost completely energeti- 
cally stabilized by the lower potential energy of the water 
molecules in the hydration shell. Consequently, the in- 
creasingly stably water-water hydrogen bonds close to a 
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hydrophobic particle, leading the positive heat capacity 
of solvation, which is a signature for hydrophobic hydra- 
tion is the key to the observed behavior. The strong 
energy and volume differences indicates a steep posi- 
tive slope of the corresponding transition line of about 
100 bar K^^. The observed stability line for the collapsed 
hydrophobic chain shows strong similarity with the lower 
temperature side of the stability digram of proteins in 



aqueous solution. 
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